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Performance of Thermal Mass Flow Meters in a 
Variable Gravitational Environment 


John E. Brooker and Gary A. Ruff 
National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 

The performance of five thermal mass flow meters, MKS Instruments 179 A and 258C, Unit Instruments 
UFM-8100, Sierra Instruments 830L, and Hastings Instruments HFM-200, were tested on the KC-135 
Reduced Gravity Aircraft in orthogonal, coparallel, and counterparallel orientations relative to gravity. 
Data was taken throughout the parabolic trajectory where the g-level varied from 0.01 to 1.8 times normal 
gravity. Each meter was calibrated in normal gravity in the orthogonal position prior to flight followed by 
ground testing at seven different flow conditions to establish a baseline operation. During the tests, the 
actual flow rate was measured independently using choked-flow orifices. Gravitational acceleration and 
attitude had a unique effect on the performance of each meter. All meters operated within acceptable limits 
at all gravity levels in the calibrated orthogonal position. However, when operated in other orientations, 
the deviations from the reference flow became substantial for several of the flow meters. Data analysis 
indicated that the greatest source of error was the effect of orientation, followed by the gravity level. This 
work emphasized that when operating thermal flow meters in a variable gravity environment, it is critical 
to orient the meter in the same direction relative to gravity in which it was calibrated. Unfortunately, there 
was no test in normal gravity that could predict the performance of a meter in reduced gravity. When 
operating in reduced gravity, all meters indicated within ± 5 % of the full scale reading at all flow 
conditions and orientations. 


INTRODUCTION 

Thermal mass flow meters (MFM) and controllers are 
commonly used in industry, laboratories, and recently in 
space flight experiments such as Combustion Module- 1, 
which first flew in 1997 aboard the space shuttle 
Columbia. These applications depend upon accurate 
fluid flow measurements to provide necessary 
information for process control, fluid delivery, and 
subsequent data analysis. In a terrestrial environment, a 
properly installed and calibrated device should be 
adequate for most applications. However, in an 
environment where the gravitational level is reduced or 
variable, special precautions should be taken to ensure 
proper delivery of gases. (Why these flow meters may 
be particularly susceptible to gravitational effects will 
be identified in a later section when discussing the 
theory of operation.) The purpose of this research is to 
determine the effect of gravitational acceleration on the 
performance of five thermal MFMs from four 
manufacturers 1 as identified in Table 1. These meters 
were selected for evaluation because they are 
representative of those used in the Microgravity Science 
Division at NASA GRC, not because they are 
representative of all commercially-available flow 
meters. No endorsement for use of these meters is 
implied. This study is intended to give the reader 


1 Only one mass flow meter from each lot was tested, therefore the 
results do not represent the entire meter population. 


information on the application of thermal MFMs and 
provide those that have used these types of instruments 
an example of potential systematic errors that may exist 
through off-nominal use. 

Table 1. Mass Flow Meters and Manufacturers 



THEORY OF OPERATION 
HFM200 

As reported by the manufacturer [1], the Hastings HFM- 
200 series instrument operates on a thermal electric 
principle whereby flow through a metallic capillary 
tube, which is a fixed percentage of the total flow 
through the instrument, is heated uniformly by a 
resistance winding attached to the midpoint of the 
capillary. Two thermocouples are welded at equal 
distances from the midpoint and develop equal outputs 
at zero flow. When flow occurs through the tubing, heat 
is transferred from the tube to the gas on the inlet side, 
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and from the gas back to the tube on the outlet side 
creating an asymmetrical temperature distribution. The 
thermocouples sense the change in capillary tube 
temperature and produce approximately a 1 millivolt 
hill scale output signal proportional to that change. This 
signal is amplified by the meter circuitry to provide a 
0-5 VDC output. For a constant power input, the 
differential thermocouple output is a function of the 
mass flow rate and the heat capacity of the gas. Since 
the heat capacity of many gases is relatively constant 
over wide ranges of temperature and pressure, the flow 
meter may be calibrated directly in mass units for those 
gases. The HFM200 sensor measures approximately 
10 seem full scale flow. Measurement of flow rates 
higher than the 10 seem full scale is achieved by 
dividing the flow with a fixed ratio shunting 
arrangement. This is accomplished by placing the 
measuring capillary tube in parallel with a laminar flow 
element that generates a pressure drop proportional to 
the mass flow rate. Therefore, the sensor only needs to 
heat the gas passing through the capillary tube. 

MKS179A, MKS258C, UFM8100, SI830L 

The theory of operation for the MKS179A, MKS258C, 
UFM8100, and SI830L is basically the same as 
described above with several exceptions [2]. At the 
midpoint of the capillary tube, two wire coils are 
wrapped side by side. The windings serve as both 
heaters and temperature sensors. The MFM electronics 
provide a constant current to the coils which are heated 
by the resistance of the wire. Since the resistance of the 
coils varies with temperature, the coils function as 
resistance temperature detectors which respond to 
changes in the gas temperature. Each coil is part of a 
resistive leg in a Wheatstone bridge circuit. When there 
is no gas flow, heat from the coils generates a uniform 
temperature gradient about the midpoint of the tube. 
When gas flows through the sensor, the upstream coil 
will cool as heat is transferred to the gas. At the 
downstream coil, the gas temperature is higher so the 
downstream coil cools only slightly. Assuming all other 
heat losses for the coils are the same, the temperature 
difference between the coils can be linearly correlated to 
the mass flow rate. The electronics in the MFM convert 
the temperature difference to a 0-5 VDC output. While 
the MKS179A, MKS258C, UFM8100, and SI830L 
instruments operate on the same principle, there are 
differences in the construction and circuitry that could 
affect their behavior in a variable gravity environment. 

Gravitational Effects 

Any instrument in which heat transfer is essential to 
obtain a measurement could be affected by variations in 
the magnitude or direction of the gravity vector. Tison 
[3] performed a detailed investigation of the 


performance of thermal flow meters and quantified the 
errors in the zero and span readings associated with the 
use of gas correction factors, meter orientation, and 
variations of upstream pressure and temperature. In this 
work, he found that the orientation of the meter in 
normal gravity changed the zero indication by less than 
0.44% of full scale and the effect on the span was 
negligible. Of course, all of tests were conducted in 
normal gravity so only the direction of the gravity 
vector was varied. Kondo et al. [4] and Weinheimer and 
Ridley [5] have investigated the effect of ambient 
pressure on thermal mall flow controllers and measured 
flow errors of factors of 2 to 3 when the pressure is 
reduced to a few millibars. These errors became 
significant only when the ambient pressure was reduced 
below about 50 mbar. Furthermore, Weinhiemer and 
Ridley [5] found that the sign of the error was reversed 
for flow meters from different manufacturers and 
depended on their specific design and operating 
characteristics. They concluded that variations in the 
heat transport involving the air near the sensing 
thermistors played a role in these errors. 

Given that thermal mass flow meters respond to any 
external stimuli that alter the heat transfer within the 
device, it is likely that variations in both the magnitude 
and direction of the gravity vector could also produce 
errors in the flow measurement. Because of the small 
diameter of the capillary tube, buoyancy plays a 
negligible role on the fluid flowing through the tube. [3] 
However, all the instruments operate by heating a 
portion of the capillary tube. A change in the natural 
convection from the outside of the capillary tube by a 
change in the gravity level or orientation of the meter 
would be interpreted by the instrument as a change in 
the flow, resulting in an inaccurate measurement. This 
effect may be alleviated by judiciously insulating the 
tube but the specific type and placement of the 
insulation could produce varying performance, even 
between otherwise identical instruments. Other 
structural differences, such as how the heating wires are 
bonded to the capillary, could also introduce an effect of 
either gravity or orientation. The purpose of this study is 
to quantify these effects. A detailed evaluation of the 
design of each MFM to determine why they perform the 
way they do was beyond the scope of this investigation. 
As previously stated, no endorsement of any flow meter 
or manufacturer is implied by NASA or the authors. 

HARDWARE DESCRIPTION 

The experiment package consisted of a test platform, a 
gas delivery system, and a data acquisition and control 
(DAC) system. On the test platform, all five mass flow 
meters were oriented in the same direction and secured 
to a plate. A critical flow orifice was installed upstream 
of each meter. The orifice provided a non-thermal flow 
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measurement with which the results could be compared. 
The platform could be positioned in different flow 
orientations relative to the gravity vector as depicted in 
Figure 1. A three-axis accelerometer was mounted on 
the plate to provide measurements of local acceleration. 


Co pa ra II el Counter para! lei O rt hog on al 



Figure 1. MFM position relative to gravity vector 


The gas delivery system contained a gas bottle (Air 
Products UHP Zero Grade Air) and various components 
to provide the test devices with a smooth and 
continuous supply of air. A flow schematic is shown in 
Appendix A. Data acquisition and control was 
performed by a PC-based measurement system as 
described below. 

The transducers used for this experiment included five 
MFMs, two pressure sensors (ambient and gas), two 
K-type thermocouples (ambient and gas), and a three- 
axis accelerometer producing a total of 12 analog 
signals. A portable, modular signal conditioning system 
(National Instruments SCC-2345) containing a thermo- 
couple input module and five lowpass filter modules 
was used in conjunction with a multifunction I/O card 
(National Instruments DAQCard-AI-16XE-50) with 
16-bit resolution and 20 kHz sampling rate for signal 
measurement and conditioning. National Instruments 
Lab VIEW Version 5 was used as the application 
software on a Dell D266XT notebook computer with 
Microsoft NT 4.0 as the operating system. Data was 
collected at 100 Hz sampling rate per channel. 

PROCEDURE 

Prior to testing, each instrument used in this experiment 
was calibrated by the NASA GRC Calibration. 
Appendix B summarizes the first-order linear regression 
models for each instrument used in this experiment. All 
mass flow meters were calibrated only in the orthogonal 
attitude in order to evaluate its performance in an off- 
nominal orientation, i.e., using the meter in an 
orientation different than the calibration procedure. 
Orifice flow rates were calculated from a first-order 
linear regression model also supplied by the Calibration 
Lab. Ground-based experiments were performed prior 
to flight to a baseline performance at normal gravity. 
Data was collected throughout the entire flight 


maneuver encompassing a variable gravitational 
environment between 0.01 and 1.8 times normal gravity. 
The analysis was performed on data sampled at one 
second increments from a ten second time interval in 
each of the three gravity regimes (normal, reduced, and 
elevated). 

DATA REDUCTION 

The data set consisted of three independent variables: 
inlet gas pressure (7 levels), gravity level (3 levels), and 
orientation (3 levels), with one dependent variable, flow 
rate. MFM performance was measured by calculating 
the difference between orifice and meter flow rate using 
the equation 



where D m f m is the meter departure function and F or f and 
F m f m are the flow rate measured by the orifice and 
thermal flow meter, respectively. Therefore, a departure 
function could be determined as a function of inlet gas 
pressure for each meter at any gravity level and 
orientation. These results will be discussed later after 
the experimental uncertainties associated with this 
experiment have been addressed. 

UNCERTAINTY ANALYSIS 

The performance of the thermal flow meters in a 
variable-g environment can only be quantified relative 
to the uncertainty existing in the measurement. The 
analysis of uncertainties in experimental measurements 
and calculations applied in this study was consistent 
with the methods outlined in the ANSI/ASME Standard 
PTC 19.1 [6] and conducted according to the procedures 
formulated in Coleman and Steele [7]. A brief summary 
of the method applied in this work is provided below. 

Propagation of Systematic Uncertainties into 
Experimentai Resuit 

Each value of an individual variable that is measured 
contains elemental error sources in the form of 
systematic and random errors. These errors propagate 
into the data reduction equation as shown in Fig. 2, and 
yield the systematic and random uncertainties of the 
experimental result. The general data reduction equation 
is given by 

r=r(X l9 X 29 ... 9 Xj) (2) 

where r is the experimental result and Xj is a measured 
variable. The systematic uncertainties from elemental 
error sources must be estimated and combined to form 
the estimate of the systematic uncertainty for each 
measured variable. The root-sum-square combination of 
M elemental errors sources for measured variable J is 
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where Bj is the systematic uncertainty and ( Bj) k is the 
uncertainty from an elemental error source. These 
systematic uncertainties propagate into the data 
reduction equation as described by 
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and the standard deviation of the average result from the 
M individual tests is 
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( 8 ) 


From these equations, the random uncertainty of the 
mean result for M > 1 0 is calculated using 


where B r is the systematic uncertainty for the result, B t 
is the systematic uncertainty of measured variable X h 
and B ik is the covariance estimator accounting for 
correlated systematic uncertainties between measured 
variables (variables that share the same elemental error 
sources). A satisfactory approximation for the correlated 
systematic uncertainty is 

L 

B * = 'L( B i)j! t k) a (5) 

a=l 


where L is the number of elemental systematic error 
sources that are common for measurements of variables 
X t and X k . 



Elemental error sources 

Individual measurement 
systems 

Measurements of 
individual variables 


Data reduction equation 


Experimental result 


P~r=tS~ r (9) 

where the value of t is taken from a t-distribution table 
at the 95% confidence interval. The systematic 
uncertainty associated with the mean result from Eq. (6) 
is the same as that found for a single test since 
systematic errors are fixed and not affected by 
averaging the results of multiple tests. Therefore, the 
total uncertainty in the mean result at the 95% 
confidence level is expressed as 

Uf = B; + P-- (10) 


First-Order Linear Regression Uncertainty 

When experimental results are represented by a 
regression model, the uncertainty associated with the 
original data will propagate into the values that are 
predicted by the model. For this analysis, a first-order 
linear model was utilized. By the method of least 
squares, the coefficients of the first-order model 
generated from N data pairs of (X h f) are 


N N N N 

Zv'Z^-ZvZv^ 


n — i=l 

Llr\ — 


i=l 


i= 1 i= 1 


N^Xf- XV 


i= 1 


XI 


\i = 1 


(ii) 


Figure 2. Propagation of errors into experimental 
results [7] 


Random Uncertainty of Experimental Result: 
Multiple Tests Method 

Calculating the random uncertainty associated with a 
data reduction equation is best accomplished from direct 
determination using multiple tests. For a single test 
replicated M times, the mean value of the result is 


r = 



( 6 ) 


where r k is the measured result from a single test. The 
standard deviation of the sample population of the M 
individual tests is given by 


N N N 

^Zv-^-Z^Z Y t 

i=l i= 1 i=l 

N (N V 

a^Z x f- Z v 

i= 1 \i=\ J 


( 12 ) 


The expression for the uncertainty in Y determined by 
the first-order regression model for a new value of X is 
found by substituting Eqs. (11) and (12) into 

Y ( X new) = a 0 + a l X new (l 3 ) 

and applying the propagation equations. As a result, the 
general form of the expression for the uncertainty in the 
predicted value of Y is 
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MFM at all three orientations and gravity levels. 
Individual plots of the results, shown in Fig. 3, will be 
discussed in the following sections. The systematic 
elemental error sources are listed in Table D.l of 
Appendix D. The values were either estimated or 
supplied by the manufacturer. With this information, 
Eq. (4) was used to calculate the propagated systematic 
uncertainties. Random uncertainties of the measured 
variables were found by the method of multiple tests 
and calculated using Eq. (9). The results are shown in 
Table D.2 of Appendix D. The total uncertainty of the 
measured variables and values determined from first- 
order linear regression models were calculated using 
Eqs. (10) and (14), respectively. The total uncertainty 
for each measurement is shown as error bars on the 
departure functions in Fig. 3. An analysis of variance 
showing the relative importance for each variation 
source was calculated from Eq. (15) and will be 
discussed later. The combined results of all calculations 
are presented in Appendix E. These results will be 
presented and discussed in the following text. 


The first seven terms account for uncertainties from the 
original (X h Y,) data pairs. The eighth and ninth terms are 
from random and systematic uncertainties associated 
with the new X value. The tenth term accounts for 
correlated systematic errors between the new X variable 
and the original X t variables. The last term accounts for 
correlated systematic errors between the new X variable 
and the original Y t variables. 

ANALYSIS OF VARIANCE 

The analysis of variance (ANOVA) was used to uncover 
the main and interaction effects of the independent 
variables (inlet gas pressure, gravity ratio, and attitude) 
on the dependent variable (flow rate). A main effect is 
the direct effect of an independent variable on the 
dependent variable. An interaction effect is the joint 
effect of two or more independent variables on the 
dependent variable. For this experiment, a three-factor 
ANOVA with replication was performed [8]. The 
relative importance of a variation source was calculated 
by using 


— s v ' 
/= l 

where MS t is the mean square of a variation source. 

RESULTS 

A typical response of the mass flow meters to 
gravitational acceleration at 50% of full scale flow in 
the three orientations is shown in Appendix C. Using 
Eq. (1), the departure function was calculated for each 


The criteria for acceptable performance used in this 
study was that the departure from orifice flow remain 
within ±1% of full scale. This was consistent with the 
manufacturer’s specification of accuracy, which 
includes the effects of non-linearity, hysteresis, and 
non-repeatability. 

In the following sections, the performance of each flow 
meter will first be discussed individually. Recall that a 
positive departure function indicates that the MFM 
reading was less than the orifice flow meter; a negative 
departure function means the MFM was reading higher 
than the orifice (see Eq. (1)). Also, because all MFM 
were calibrated in the orthogonal configuration in 
normal gravity, the best performance would be expected 
at this orientation and gravity level. In Fig. 3, note that 
the scale on the vertical axis has been adjusted for each 
meter to maintain similar graphical resolution. 

After discussing the performance of each meter, the 
results of the analysis of variance will be presented, 
followed by a discussion of the overall performance of 
the mass flow meters. 

MKS179A 

In normal gravity, the departure function, i.e., the 
difference between the flow rate measured by the orifice 
and that measured by the MKS179A, differed by less 
than ±1% of full scale when the flow meter was in the 
orthogonal or counterparallel orientations. However, in 
the coparallel configuration, the departure is greater 
than 1% of full scale and positive, indicating that the 
MFM was reading less than the orifice flow meter. If 
convective heat transfer plays a role in producing these 
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variations in 1-g, variations with gravity level would 
also be expected. This was indeed observed, as shown 
in Fig. 3. However, at both elevated and reduced 
gravity, the departure function was within ±1% of full 
scale when the MFM was in the orthogonal (calibrated) 
orientation. When oriented coparallel, the MKS179A 
indicated less than the orifice and when counterparallel, 
it indicated greater than the orifice. In elevated gravity, 
the departure was greater than ±1% of full scale for all 
but the lowest flow rate; in reduced gravity, the 
departure increased with increasing flow rate but was 
greater than ±1% of full scale only for the highest flow 
rates. 

MKS258C 

This meter is also manufactured by MKS so its behavior 
could be expected to be similar to the MKS 179 A. 
However, as shown in Fig. 3, there are significant 
differences. In the orthogonal orientation, the departure 
from the orifice reading is within ±1% of full scale for 
normal and reduced gravity level. For elevated gravity 
levels, it is within ±1% of full scale for all but the 
highest flow rates. Co- and counter-parallel orientations 
at elevated and normal gravity showed large departures 
from the orifice reading with the MKS258C reading 
higher than the orifice in the counterparallel orientation 
and lower than the orifice in the coparallel orientation. 
The large deviation in flow rate with orientation even in 
normal gravity emphasizes the importance of using a 
flow meter in the orientation in which it was calibrated. 
In spite of this, the departure functions for all 
orientations are essentially within ±1% of full scale in 
reduced gravity. This behavior would be expected in the 
absence of gravity if changes in natural convection was 
the only cause of the deviations observed in non-zero 
gravity. 

UFM8100 

The UFM8 1 00 demonstrated unique behavior in several 
respects. First, even in the orthogonal (calibrated) 
orientation at normal gravity, the departure function is 
greater than ±1% for almost all flow rates. This is 
indicative of a calibration problem with either the mass 
flow meter or the reference orifice. Unfortunately, the 
experimental apparatus was disassembled before this 
data was analyzed so the cause of this behavior could 
not be verified. However, this appears to be a systematic 
effect that could be corrected upon re-calibration. 
Figure 3 shows the UFM8100 meter demonstrated 
behavior within about 1.5% of foil scale of the 
orthogonal departure for all orientations. Assuming that 
conclusions can be drawn based on the effect of attitude 
and g-level, this indicates that the output of this meter 
was nearly independent of orientation and gravity level. 


As with the MKS258C, there was little influence of 
orientation in reduced gravity. 

SI830L 

For the SI830L, note that the vertical scale is 
considerably larger than the other instruments, making 
the ±1% of full scale band much narrower. In an 
orthogonal configuration in normal gravity, i.e., the 
configuration in which it was calibrated, the departure 
function for the SI830L is within ±1%FS. However, if 
operated in an off-nominal orientation in normal 
gravity, the departure function have absolute values of 
nearly 10% of full scale. In the coparallel orientation, 
the SI830L reads lower than the orifice while in the 
counterparallel orientation, it reads higher, similar to the 
behavior of the MKS 179 A and MKS258C meters. The 
behavior of the SI830L at elevated gravity levels is 
similar to its behavior in normal gravity. In reduced 
gravity, there is little effect of orientation on the 
departure function although all three orientations have 
departure functions greater than ±1% of full scale with 
the SI830L reading greater than the orifice flow meter. 
Even though there is little effect of orientation on flow 
rate in reduced gravity, there appears to be an overall 
effect of reduced buoyancy on the output of the meter. 
(This can also be seen in Appendix C in the way the 
trace of the SI830L output closely mirrors that of the 
vertical accelerometer.) 

HFM200 

The behavior of the HFM200 meter is similar to that of 
the MKS 179 A meter. The difference is that at normal 
gravity, the departure functions are within 1%FS for all 
orientations at all inlet gas pressures. Since there is no 
effect of orientation on meter output at normal gravity, 
one might expect the output to also be independent of 
orientation in reduced gravity. However, as shown in 
the figure, there is a significant effect with the coparallel 
configuration having departure functions greater than 
±1% of full scale at all inlet pressures. The departure 
function for the counterparallel configuration is greater 
than ±1% of foil scale only at the highest inlet pressure. 
The reason for this behavior was not identified 
although, given the results in normal gravity, it would 
not appear that it was solely caused by a difference in 
natural convection. 

ANOVA Results 

An ANOVA (analysis of variations) was performed on 
these data to quantify the relative importance of changes 
in orientation, gravity level, inlet gas pressure, and the 
interactions between these variables on the values of the 
departure functions. These results are shown in Table 2. 
For all meters except the UFM8100, the meter 
orientation had the largest effect on the departure 
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function independent of gravity level. This emphasizes 
the importance of using a mass flow meter in the 
orientation at which it was calibrated. For these meters, 
the gravity level was the second highest source of 
variation but this was generally an order of magnitude 
less significant than the orientation. 


Table 2. Three-Factor ANOVA for Departure Function 


Meter 

A 

% Relative Importance 
B C AB AC BC 

ABC 

MKS179A 

86.9 

8.92 

0.25 

2.95 

0.84 

0.06 

0.06 

MKS258C 

78.2 

2.60 

1.62 

17.1 

0.21 

0.20 

0.09 

UFM8100 

13.7 

23.9 

58.7 

2.15 

1.05 

0.23 

0.20 

SI830L 

79.6 

1.81 

0.12 

18.4 

0.02 

0.02 

0.02 

HFM200 

79.2 

2.40 

0.16 

16.6 

1.21 

0.06 

0.33 


A - Orientation 
B - Gravity Level 
C - Inlet Gas Pressure 


The inlet gas pressure had the largest contribution to the 
departure function for the UFM8100 meter which is 
indicative of a unit that is out of calibration, as 
previously discussed. Assuming this is a known source 
of error that could be eliminated upon re-calibration, the 
next largest contribution to the error for this meter was 
the gravity level, not the orientation as for the other 
meters. This is unique among the meters tested and 
indicates that the UFM8100 had the least sensitivity to 
orientation at all gravity levels. For all meters, the effect 
of joint interactions between inlet gas pressure, 
orientation, and gravity level were much smaller than 
the major sources of variation and can be neglected. 

DISCUSSION 

An ideal flow meter would have an output that was 
independent of orientation and gravity level. In terms of 
the departure functions shown in Fig. 3, all the curves 
for this ideal meter would be coincident and lie within 
the ±1% FS band. Similar to previous investigations, 
this study showed flow measurement errors that were 
dependent on orientation of the meter as well as the 
magnitude and direction of the gravity vector. Tison [3] 
had found that the orientation of the meter in normal 
gravity changed the zero indication by less than 0.44% 
of full scale and the effect on the span was negligible. 
The results from the HFM200 and UFM8100 are 
consistent with this finding. However, the variation in 
the zero level with orientation in normal gravity ranged 
from about ±2.5% of full scale for the MKS179A and 
MKS258C meters and up to ±10% of full scale for the 
SI830. The MKS179A meter was the only one for 
which a significant error in the span was observed. The 
current results are consistent with the those of the 


previous investigations [3-5] in which variations were 
observed between meters from different manufacturers. 
Also, since Tison [3] evaluated meters having 
considerably lower flow ranges than the meters 
evaluated in this study, it is plausible that there could be 
significant variations in the magnitude of the errors. 

If differences in buoyancy and therefore, natural 
convection heat transfer, were the only factor causing 
the output of thermal flow meters to vary, the effect of 
orientation would become small in reduced gravity 
since the g-level is similar in any direction. This type of 
behavior was observed for the MKS258C, UFM8100, 
and SI830L instruments but not for the MKS179A or 
HFM200. Also, if a change in output was caused by 
buoyancy, the change should become greater in elevated 
gravity since buoyancy is enhanced in this environment. 
This behavior was observed for the MKS258C and 
SI830L; the other meters showed approximately the 
same behavior in elevated gravity as in normal gravity. 
Surprisingly, the UFM8100 meter showed similar 
behavior in all gravity levels with little effect of 
orientation. 

As noted above, the output of the MKS179A and 
HFM200 instruments was dependent upon orientation 
even in reduced gravity. This indicates that more 
complex heat transfer mechanisms exist between the 
internal components of these instruments. While a 
detailed evaluation of the design of each flow meter and 
its heat transfer characteristics was beyond the scope of 
this study, there are several potential explanations. First, 
the time scale of the thermal response could be longer 
than 25 seconds so the meter did not reach a thermal 
steady-state during the period of reduced gravity. 
Alternatively, the flow meter could be sensitive to the 
high-frequency, low amplitude fluctuations in g-level, 
z.e., g-jitter, that occurs during the low- gravity portions 
of parabolic flight. The former mechanism was assumed 
to be negligible because the response time of all 
instruments was observed to be fast relative to the 
duration of the low gravity period (see Appendix C). 
However, the influence of oscillations in g-level was not 
specifically addressed and could be a factor. This effect 
could be evaluated in a drop tower where the residual g- 
level is lower although maintaining instrument 
calibration and repeatability through the high g-levels 
experienced at the end of a drop would require special 
care. 

The data obtained in this study also shows that there is 
no test in normal gravity that will yield conclusive 
information about how a specific thermal flow meter 
will respond in either reduced or elevated gravity For 
four of the five flow meters, the variation in the 
departure function with inlet gas pressure at elevated 
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gravity is similar to the variation observed in normal 
gravity. Therefore, normal gravity behavior is a 
reasonable indicator of meter performance in elevated 
gravity. Unfortunately, these results show that the 
variation of departure function with inlet pressure at 
reduced gravity may be the same (MKS179A), better 
(MKS258C, UFM8100, SI830L), or worse (HFM200) 
than the behavior in normal gravity. In all cases, 


however, the departure functions at reduced gravity 
were within about ±5% of full scale. This indicates that 
all of these meters could be used in reduced gravity if a 
higher uncertainty is accepted. During operation in long 
periods of reduced gravity, it is recommended that 
measures be made to periodically check the calibration 
of thermal flow meters. 
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Figure 3. Departure from orifice flow for three flow meter orientations and gravity levels 
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Figure 3. Departure from orifice flow for three flow meter orientations and gravity levels (concluded) 




SUMMARY AND CONCLUSIONS 

The performance of five thermal mass flow meters were 
tested on the KC-135 Reduced Gravity Aircraft in 
orthogonal, coparallel, and counterparallel orientations 
relative to the g vector. Data was taken throughout the 
parabolic trajectory where the g-level varied from 0.01 
to 1.8 times normal gravity. Each meter was calibrated 
in normal gravity in the orthogonal position prior to 
flight followed by ground testing at seven different flow 
conditions to establish a baseline operation. Departure 
functions relative to an orifice flow reference was 
developed to quantify the results. A three-factor analysis 
of variance was performed to uncover the main and 
interaction effects of the inlet gas pressure, gravity 
level, and flow meter orientation. A detailed uncertainty 
analysis was also performed to quantify these effects 
relative to the systematic and random errors of the 
experiment. 

All of these flow meters performed within 
manufacturers specifications when operated in normal 
gravity in the orientation at which they were calibrated. 
This, of course, is how they were designed and is 
recommended by the manufacturer. Gravitational 
acceleration and orientation had a unique effect on the 
performance of each meter. Specific conclusions from 
this study are as follows: 

1. For only three of the five flow meters did the 
departure functions for the different orientations 
converge in reduced gravity. This implies that more 
complex heat transfer mechanisms exist between the 
internal components of these instruments than could 
be evaluated in this experiment. Explanations for this 
behavior include (1) a transient thermal response 
such that the meter did not reach steady state during 
the 25 seconds of reduced gravity on the aircraft or 

(2) sensitivity to the high-frequency, low amplitude 
fluctuations in g-level that occurs during the low- 
gravity portions of parabolic flight. 

2. For four of the five meters, the variation in the 
departure function with orientation at normal gravity 
was similar to that experienced at elevated gravity. 
Therefore, meter performance in normal gravity 
gives a reasonable indication of its performance at 
elevated gravity levels. 

3. The ANOVA analysis indicated that for four of the 
five meters, variations in the orientation had a 
greater contribution to the departure function than 
variations in gravity level. This indicates that if a 
flow meter will be used in a non-zero gravity 
environment, the meter should be mounted in the 
attitude at which it was calibrated. 

4. There is no single test that can be performed in 
normal gravity that will conclusively indicate the 
reduced gravity performance of a thermal flow 


meter. Based on this evaluation of five flow meters, 
the variation of departure function with inlet 
pressure at reduced gravity may be the same, better, 
or worse than the behavior in normal gravity. 

5. In reduced gravity, the performance of all meters 
evaluated was within about ±5% of full scale for any 
orientation, independent of its performance in 
normal or elevated gravity levels. 

6. Even though the flow meters operate on similar 
principles, their performance was quite different, 
presumably because of variations in the internal 
design and details of operation. This emphasizes the 
need for care in handling thermal flow meters after 
they have been calibrated. Small changes in 
configuration, mounting, or internal structure could 
have dramatic and unexpected effects on their 
operation and accuracy. 

Based on the experience gained during this 
experiment and the subsequent results, it is suggested 
that when operating thermal flow meters in low-gravity 
environments, the greatest accuracy can be obtained by 
the in- situ calibration of the flow meter. 
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APPENDIX A: FLOW SCHEMATIC 
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15 PSI A 
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Figure A. 1 Schematic of the flow system used in this experiment 
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APPENDIX B: CALIBRATION AND MODELS 


Table B. 1 Calibration Reports with First-Order Linear Regression Models 


Device # 

Report # 

ao 

Value 

Units 

Value 

ai 

Units 

13 

1000046378 

-0.005 

psia 

19.995 

psia/VDC 

15 

421935 

-0.01107 

slma 

0.01047 

slma/psia 

16 

421934 

-0.00993 

slma 

0.01038 

slma/psia 

17 

421933 

-0.01081 

slma 

0.01035 

slma/psia 

18 

421932B 

-0.01557 

slma 

0.01057 

slma/psia 

19 

421 931 C 

-0.01025 

slma 

0.01080 

slma/psia 

20 

1000048273 

5.559E-03 

slma 

0.1987 

slma/VDC 

21 

1000048272 

1.683E-03 

slma 

0.1998 

slma/VDC 

22 

1000048271 

9.158E-04 

slma 

0.1997 

slma/VDC 

23 

1000048270 

-2.132E-03 

slma 

0.1984 

slma/VDC 

24 

1000048269 

6.306E-03 

slma 

0.1979 

slma/VDC 

25 

1000038626 

-2.338 

psia 

3.018 

psia/VDC 


Form of Equation: Value = 3,1 (Output) + a 0 
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APPENDIX C: RESPONSE OF THE MASS FLOW METERS 
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Figure C.l Response of the MFM gravitational acceleration and attitude (50 %FS) 
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APPENDIX D: ERROR SOURCES AND UNCERTAINTIES 


Table D.l Systematic Elemental Error Sources 


Variable 

Source 

Notes 

F mfm 

Calibration 

DAC 

0.2% of reading or 0.02% FS 
(1 slma), whichever is greater 
Manufacturer-supplied equation 

Pg 

Calibration 

DAC 

0.03% FS (100 psia) 
Manufacturer-supplied equation 

Pa 

Calibration 

DAC 

0.59% FS (15 psia) 
Manufacturer-supplied equation 

Tg 

Calibration 

DAC 

1.0 K (manufacturer) 
1.0 K (manufacturer) 

T a 

Calibration 

DAC 

1.0 K (manufacturer) 
1.0 K (manufacturer) 


g Calibration 2% FS (2 g 0 ) 


DAC Manufacturer-supplied equation 


Table D. 2 Random Uncertainties 


Variable 

Ground 

(normal) 

Aircraft 

(reduced) 

Aircraft 

(elevated) 

Fmksi79a (slma) 

0.000 

0.001 

0.001 

Fmks258c (slma) 

0.000 

0.001 

0.001 

Fufmsioo (slma) 

0.000 

0.002 

0.001 

F si83ol (slma) 

0.000 

0.003 

0.004 

Fhfm2oo (slma) 

0.000 

0.001 

0.001 

Pg(psia) 

0.04 

0.08 

0.07 

P a (psia) 

0.00 

0.05 

0.04 

T g (°C) 

0.0 

0.5 

0.5 

Ta(°C) 

0.0 

0.7 

1.0 

g z /go 

0.00 

0.01 

0.02 
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APPENDIX E: TABULATED RESULTS 


The tables in this appendix contains the raw and reduced data for all five mass flow meters, gravity levels, and flow 
meter attitudes. In all tables, parentheses around numbers indicates a negative value. 
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Mean 0.619 0.611 0.603 0.733 0.598 56.92 12.23 12.4 15.0 (0.04) 1.77 (0.09) 0.604 0.586 0.578 0.581 0.585 (1.5) (2.5) (2.5) (15.2) (1.3) 

± 0.001 0.001 0.001 0.002 0.001 0.05 0.09 1.4 1.4 0.04 0.04 0.04 0.001 0.001 0.001 0.001 0.001 0.2 0.2 0.2 0.2 0.2 


APPENDIX F: SYMBOLS 


B 

i 

Bj 

B r 


DAC 

D 

F 

FS 

g 

go 

g z 

Hz 

J 

L 

M 

MFM, mfm 

MSi 

N 

orf 


P? 

r 

f 

Rli 

s 

r 

Sj 

seem 

slma 

t 

T 

a 

T 


coefficients of linear regression model 

estimate of systematic uncertainty at a 95% confidence limit for i th measured variable 
systematic uncertainty for measured variable J 
systematic uncertainty of result 

covariance estimator for L correlated (common) elemental systematic error sources for 

measurements of variables X t and X k 

data acquisition and control 

departure function 

mass flow rate 

full scale 

gravitational acceleration vector 
vertical gravitational acceleration at surface 
vertical gravitational acceleration 
hertz 

total number of measured variables 

total number of elemental systematic error sources common for measurements of variables X t and X k 
total number of elemental error sources; total number of replications 
mass flow meter 

mean square of a variation source 

number of data pairs 

orifice 

ambient pressure 
gas pressure 

random uncertainty of mean result 

experimental result determined from J measured variables 

mean value of result from replicating a single test M times 

relative importance of a variation source 

standard deviation of sample population of M individual tests 

standard deviation of mean result from M tests 

standard cubic meters per minute air 

standard liters per minute air 

t-distribution 

ambient temperature 

gas temperature 
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VDC volt direct current 

U- total uncertainty of mean result 

U Y total uncertainty of y- value determine from the first-order linear regression model 

X. i th measured value of variable; x- value of i th data pair 

X new value of X used in 1 st -order linear regression model 

Y(X new ) l st -order linear regression model 

Y. y- value of i th data pair 
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